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Abstract 

In this paper, we prove that if the matrix of the linear system is sym- 
metric, the Cholesky decomposition can be obtained from the Gauss elim- 
ination method without pivoting, without proving that the matrix of the 
system is positive definite. 
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1 Introduction 

The direct methods for solving linear system , Gauss elimination method, LU 
decomposition and Cholesky method are well known. We agreed with some 
authors [T] , [2] that the LU decomposition and the Cholesky method are helpful 
for solving many linear systems of the same matrix, whose difference is only the 
constants at the right hand side. 

The Gauss elimination method with or without pivoting can lead us to the 
LU decomposition. The Gauss elimination method with pivoting applied to 
linear system with symmetric matrix can't lead us to the Cholesky decomposi- 
tion because pivoting alters the symmetry. However, sometimes wc don't need 
pivoting. 

Since the Gauss elimination can lead us to the LU decomposition, we think 
that it is better to solve at first one of these linear systems by the Gauss elim- 
ination and the others by LU decomposition. So, we will be allowed to solve 
at first one of these linear systems by the Gauss eliminaton and the others by 
Cholesky method, in the case where the matrix is symmetric positive definite, 
if Gauss elimination can lead us to the Cholesky method. 

In this paper, we will sec that there exist relation between Gauss elimination 
without pivoting and the Cholesky method. That is Gauss elimination without 
pivoting can lead us to Cholesky decomposition. We will prove this relation by 
using the LU decomposition. 

We organize the paper as the following. In the second section we represent 
the LU decomposition by using the Gauss elimination method. In the third 



section we will prove, with the help of the LU decomposition we can have the 
Cholesky decomposition from the Gauss elimination method. At the last section 
an example will be tread for clarify the method. 

2 Gauss elimination 

Let us consider the following system of linear n equations and n unknowns 
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can be written under matricial form 
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If an ^ 0, Gauss elimination in the first column is written 
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which may matricially be written 

AiX = Bt (2) 



and can be obtained in multiplying ([T|) by the lower triangular matrix[3], 0] 
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A 1 = GiA and B 1 = G\B. 



The Gauss elimination at the second column can be obtained by multiplying 
to the relation @ the lower triangular matrix 
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The matricial equation @ becomes 



A 2 X = B 2 



with A 2 = G 2 GiA and B 2 = G 2 G\B. 
In continuing so, we have finally 

UX = B' 
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with U = G n _i . . . G2G1A an upper triangular matrix and B' = G n —i . . . G 2 G\B. 
We can check easily that the inverse of G; is 
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So, G„_iG„_2 • • • G2G1 is an invertible matrix and the equation ([3]) becomes 

LUX = B 

with L = G^Gz 1 ■ ■ ■ G-\G-\ and we have the decomposition of A as product 
of lower triangular matrix by an upper triangular matrix 

A = LU 
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3 Cholesky method 

Now suppose A is a symmetric positive definite matrix. Then, the Cholesky 
method consists to decompose A as the product 

A = G T G 



with G is a upper triangular matrix and G T his transpose. 
Let us at first generalize this decomposition. 

Definition 1 Let A = (a.ij) 1<i j <n a complex symmetric matrix. Let us call 
Gauss matrix of A the upper triangular matrix U(A), obtained after transform- 
ing A by the Gauss elimination above. 



Proposition 2 Let A = (aij) 1<{ j <n a complex symmetric matrix, such that 
det(A) ^ 0, 
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the Gauss matrix of A. Then A can be decomposed as the product A — G T G 
with 
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with y/uu d root squared of the complex number ua. 
Proof. 
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with 



D = 



V - 7^J 

Let G = DC/ (A). Since A symmetric, hence L(A)D^ 1 = G T 
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Example 3 Consider the following two systems of linear equations 
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whose difference is only the right hand side and their matrix is symmetric. So 
let us solve the first one by the Gauss elimination method, and the second one 
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by the Cholesky method or LU method. 
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-2, X2 = 1, X\ — 1 + 1 



_B and solved in writing 
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